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INTRODUCTION 

The Van-der-Waals-type range dependence of the nuclear forces has provided a 
good part of the motivation for fragmentation studies. The predictions of a liquid- 
gas phase transition in nuclear matter [| have raised the hope that signals of it 
may be identified in reactions of finite nuclei. The observation of multifragmentation 
|| [|, ^] represented a major breakthrough in this direction, as it indicated that 
the created nuclear systems may actually pass through states of high temperature 
and low density during the later stages of the reaction. Multifragmentation has been 
predicted [f| [7| to be the dominant decay mode under such conditions that are similar 
to those expected for the coexistence region in the nuclear-matter phase diagram. 

Many features of multifragmentation are well reproduced by the statistical mul- 
tifragmentation models || ^, |j. They predominantly include the fragment distri- 
butions and correlations describing the populated partition space. But also kinetic 
observables such as the energies or velocity correlations of the produced fragments 
have been found to agree with the statistical predictions in certain cases || [10| . The 
essential assumption on which these models are based is that of a single equilibrated 
breakup state at the end of the dynamical evolution and of a statistical population 
of the corresponding phase space. To the extent that it is realized in nature, this 
scenario offers the possibility to map out the nuclear phase diagram, even though for 
finite nuclei, by sampling the thermodynamic equilibrium conditions associated with 
the multi-fragment breakup channels. 

Considerable progress has been made with this program during recent years. 

Multifragmentation has been studied over a wide range of different classes of high- 
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energy reactions and breakup temperatures have been deduced from the measured 
data. The correlation of the temperature with the excitation energy, often referred to 
as the caloric curve of nuclei and first reported for spectator decays following 197 Au 
+ 197 Au reactions at 600 MeV per nucleon [O], has also been derived for other cases 
HID] , |T2"| , [TB| ]. Critical evaluations of the obtained results and of the applied methods 



have followed. More work is clearly needed in order to complete the picture [H 



In this contribution, the focus will be on the decay of excited spectator nuclei 
15]. New results will be reported that were obtained in two recent experiments 
with the ALADIN spectrometer at SIS in which reactions of 197 Au on 197 Au in the 
regime of relativistic energies up to 1 GeV per nucleon were studied. In the first 
experiment, the ALADIN spectrometer was used to detect and identify the products 



of the projectile-spectator decay [16]. The Large- Area Neutron Detector (LAND) was 
used do measure coincident free neutrons emitted by the projectile source. In the 
second experiment, three multi-detector hodoscopes, consisting of a total of 216 Si- 
CsI(Tl) telescopes, and three high-resolution telescopes were positioned at backward 
angles to measure the yields and correlations of isotopically resolved light fragments 



of the target-spectator decay JTTfl- From these data excitation energies and masses, 



temperatures, and densities were deduced. Before proceeding to the discussion of 
these new data some of the indications for equilibration during the spectator decay 
will be briefly recalled. 

EVIDENCE FOR EQUILIBRATION 

The universal features of the spectator decay, as apparent in the observed Z bound 
scaling of the measured charge correlations, were the first and perhaps most striking 
indications for equilibrium Hl6| . The quantity Zb oun d is defined as the sum of the 



atomic numbers Zi of all projectile fragments with Zi > 2. It represents the charge 
of the original spectator system reduced by the number of hydrogen isotopes emitted 
during its decay. 

The invariance of the fragmentation patterns, when plotted as a function of Z^und, 
suggests that the memory of the entrance channel and of the dynamics governing the 
primary interaction of the colliding nuclei is lost. This feature extends to other 
observables. The transverse-momentum widths of the fragments, as shown in Fig. 1, 
do not change with the bombarding energy, indicating that collective contributions 
to the transverse motion are small. The equilibration of the three kinetic degrees 
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Figure 1: Widths of the transverse- momentum distributions cr(pt) as a function of the 
fragment atomic number Z for the reactions 197 Au on 197 Au at E/A = 600, 800, and 1000 
MeV and for 20 < Z\, oun ^ < 60. The line is proportional to yZ (from Ref. pp . 



of freedom in the moving frame of the projectile spectator was confirmed by the 
analysis of the measured velocity spectra [T(|. The square-root dependence on the 
atomic number Z implies kinetic energies nearly independent of Z and hence of the 
mass. 

The success of the statistical multifragmentation model in describing the observed 
population of the partition space may be seen as a further argument for equilibration. 
Here the main task consists of finding an appropriate ensemble of excited nuclei to be 
subjected to the multi-fragment decay according to the model prescription. Starting 
from the entrance channel may not necessarily provide sufficiently realistic ensembles, 
even though a good description of the fragment correlations was obtained with the 
quantum-molecular-dynamics model coupled to the statistical multifragmentation 
model ]nj. An alternative method consists of using empirically derived ensembles. 
Near perfect descriptions of the measured correlations, including their dispersions 
around the mean behaviour, can be achieved || [L7|. The mathematical procedure 
of backtracing allows for studying the uniqueness of the obtained solutions and their 



sensitivities to the observables that were used to generate it flT) 



The ensemble derived empirically for the reaction 197 Au on 197 Au at 1000 MeV 
per nucleon is shown in Fig. 2. Its capability of reproducing the measured mean 
multiplicity of intermediate-mass fragments and the mean charge asymmetry of the 
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Figure 2: Excitation energy E x /A as a function of the mass Aq for the ensemble of excited 
spectator nuclei used as input for the calculations with the statistical multifragmentation 
model. The area of the squares is proportional to the intensity (from Ref. fl7|| ). 



two heaviest fragments is illustrated in Fig. 3 where the dashed and dotted lines 
show the model results for E x /A chosen 15% above and below the adopted values. 
In the region Zbound > 30, the mean excitation energy of the ensemble of spectator 
nuclei was found to be well constrained by the mean fragment multiplicity alone. 
At Z bound ps 30 and below, the charge asymmetry was a necessary second constraint 
while, at the lowest values of Z bound , neither the multiplicity nor the asymmetry 
provided rigid constraints on the excitation energy. 

The spectator source, well localized in rapidity |16| and, apparently, exhibiting 
so many signs of equilibration, seems an excellent candidate for studying the nuclear 
phase diagram. Dynamical studies also support this conclusion []20], |21|j . There are 
limitations, however, which are mainly seen in the emission of nucleons and very light 
particles. Here the components from reaction stages that lead to the formation of 
the spectator system and from its subsequent breakup overlap, thereby creating diffi- 
culties for the extraction of the equilibrium properties. This is particularly apparent 
in the case of the excitation energy which will be discussed in the next section. 



EXCITATION ENERGY 



A method to determine the excitation energy from the experimental data was 
first presented by Campi et al. |22[] and applied to the earlier 197 Au + Cu data 
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Figure 3: Mean multiplicity of intermediate-mass fragments (Mjmf) (top) an d mean 
charge asymmetry (012) (bottom) as a function of Z^ oun ^ as obtained from the calcula- 
tions with the statistical multifragmentation model (open circles), in comparison to the 
experimental data for 197 Au on 197 Au at E/A = 1000 MeV (dots, from Ref. @). The 
dashed and dotted lines show the results of the calculations with excitation energies E x /A 
15% above and 15% below the adopted values, respectively. Note that the trigger threshold 
affected the data of Ref. at Zbound > 65 (from Ref. [17]). 



2~3| . It is based on the idea of calorimetry which requires a complete knowledge 
of all decay products, including their atomic numbers, masses, and kinetic energies. 
In this work, the measured abundances for Z > 2 were used and, e.g., the yields 
of hydrogen isotopes were deduced by extrapolating to Z = 1. In the same type 
of analysis with the more recent data for 197 Au + 197 Au at 600 MeV per nucleon, 
the data on neutron production measured with LAND were taken into account JTT ] . 
Since the hydrogen isotopes were not detected assumptions concerning the overall 
N/Z ratio of the spectator, the intensity ratio of protons, deuterons, and tritons, 
and the kinetic energies of hydrogen isotopes had to be made. 

The latest evaluation of the excitation energy included the measured neutron data 
for three bombarding energies and the data for hydrogen emission from the target 



spectator at 1000 MeV per nucleon [ 24[ . For the case of 600 MeV per nucleon, the 
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Figure 4: Mean excitation energy per nucleon (Eq)/{Aq) as a function of the spectator 
mass (Aq) for the reaction 197 Au on 197 Au at three bombarding energies as indicated. The 
open symbols denote the upper and lower limit of the mean excitation energy that may be 
derived from analyses with the statistical multifragmentation model (from Ref. |24j|). 



difference to the published energy values IT amounted to about 10%. More im- 
portantly, however, the deduced spectator energy was found to depend considerably 
on the bombarding energy (Fig. 4). It increases by about 30% over the range 600 
to 1000 MeV per nucleon which is in contrast to the universality observed for other 
observables (previous section) . The origin of this rise lies solely in the behavior of the 
mean kinetic energies of neutrons in the spectator frame (Fig. 5). They have a large 
effect on the deduced total excitation energy, first, because the neutron multiplicity 
is large and, second, because the hydrogen isotopes, measured only at 1000 MeV per 
nucleon, were assumed to scale in the same way as the neutrons do. 

It is obvious that this uncertainty represents a considerable problem, a memory 
of the entrance channel is inconsistent with the idea of measuring thermodynamic 
properties of an equilibrated breakup state. It is reasonable to assume that part of 
the experimentally determined energy may be due to pre-equilibrium or pre-breakup 
emission, even though the analysis of nucleon emission is restricted to the data at for- 
ward (backward) angles in the projectile (target) frame (Fig. 6). The experimental 
excitation energies are larger than the range of energies potentially consistent with 
the statistical multifragmentation model (Fig. 4), and the spectra of hydrogen iso- 
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Figure 5: Mean kinetic energy of neutrons (full symbols) in the rest frame of the projectile 
spectator as obtained in measurements including LAND for the reaction 197 Au on 197 Au 
at E/A = 600 (triangles), 800 (squares), and 1000 MeV (circles). For E/A = 1000 MeV, 
the measured mean kinetic energies of protons (open circles) are shown in comparison to 
the sum of the neutron kinetic energies and an estimated Coulomb contribution (full line, 
fromRef. @). 



topes, measured with the high-resolution telescopes, exhibit yields and slopes much 



larger than predicted by the model [|17| . The process of spectator formation involves 
secondary scatterings of fireball nucleons on spectator matter which may generate 
a pre-breakup source centered close to the spectator rapidity. Experimentally, the 
next step should consist of complementing the neutron data with equivalent data for 
proton and light-charged-particle emission at several bombarding energies. 



TEMPERATURE 

The shape of the caloric curve ||11|| , reminiscent of first-order phase transitions 
in ordinary liquids, and its similarity to predictions of microscopic statistical models 
H 0, |T7|], has initiated a widespread discussion of whether nuclear temperatures of 
this magnitude can be measured reliably (see Refs. p5| , |26|| and references given 



in these recent papers) and whether this observation may indeed be linked to a 
transition towards the vapor phase P7| , 



Here we restrict ourselves to new results obtained from the study of the target 
spectator at backward angles in the laboratory in 197 Au + 197 Au collisions at 1000 
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Figure 6: Rapidity spectra of neutrons measured with LAND for three bins in Zb oun d for 
the reaction 197 Au on 197 Au at E/A = 1000 MeV and a cut in transverse momentum of 35 
MeV/c < pt < 70 MeV/c. The thick lines represent the result of fits assuming a projectile 
and a mid-rapidity source. The individual contributions of these sources are given by the 



thin lines (from Ref. [24]). 



MeV per nucleoli. From the isotopically resolved yields of hydrogen, helium, and 



lithium isotopes breakup temperatures T He Li, ^Hepd, and T He dt were derived p9fl . The 
corrections for sequential feeding of the ground-state yields, based on calculations 
with the quantum statistical model []30|| , resulted in good qualitative agreement for 



the three temperature observables [ITT 



The corrected temperature T He Li is shown in Fig. 7. With decreasing Z bound , it 
increases from T = 4 MeV for peripheral collisions to about 10 MeV for the most 
central collisions. Within the errors, these values are in good agreement with those 
measured for projectile spectators in the same reaction at 600 MeV per nucleon. In 
both cases, the displayed data symbols represent the mean values of the range of sys- 
tematic uncertainties associated with the two different experiments while the errors 
include both statistical and systematic contributions. The projectile temperatures 
are the result of a new analysis of the original data and are somewhat higher, be- 
tween 10% and 20%, than those reported previously |Tl |. Their larger errors follow 



from a reassessment of the potential 4 He contamination of the 6 Li yield caused by Z 
misidentification. The invariance of the breakup temperature with the bombarding 
energy is consistent with the observed universality of the spectator decay. 

Calculations with the statistical multifragmentation model were performed for 
the ensemble of excited spectator nuclei shown in Fig. 2. Results have already been 
shown in Fig. 3. The solid line in Fig. 7 represents the thermodynamic temperature 
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Figure 7: Temperatures Tn e Li of the target spectator for the reaction 197 Au on 197 Au 
at E/A = 1000 MeV (dots) and of the projectile spectator at E/A = 600 MeV (open 
squares) as a function of Zb OU nd- The data symbols represent averages over bins of 10- 
units width and, for clarity, are laterally displaced by 1 unit of Zb oun d- Statistical and 
systematic contributions are included in the displayed errors. The lines are smoothed fit 
curves describing the breakup temperature Tb rea kup (full line) and the isotopic temperature 
TueLi (dashed line) calculated with the statistical multifragmentation model (from Ref. 

0)- 

T obtained in these calculations. With decreasing Zbound, & increases monotonically 
from about 5 to 9 MeV. Over a wide range of Zb OU nd it remains close to T = 6 MeV 
which reflects the plateau predicted by the statistical multifragmentation model for 
the range of excitation energies 3 MeV < E x /A < 10 MeV J/J. In model calculations 
performed for a fixed spectator mass, the plateau is associated with a strong and 
monotonic rise of the fragment multiplicities. Experimentally, due to the decrease 
of the spectator mass with increasing excitation energy (Fig. 2), the production of 
intermediate-mass fragments passes through a maximum in the corresponding range 
of Zb oun d from about 20 to 60 (Fig. 3). 

The dashed line gives the temperature T He Li obtained from the calculated isotope 
yields. Because of sequential feeding, it differs from the thermodynamic temperature, 
the uncorrected temperature TH e Li,o being somewhat lower. Here, in order to permit 
the direct comparison with the experimental data in one figure, we display Tn Li 
which has been corrected in the same way with the factor 1.2 suggested by the 
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Figure 8: Measured isotope temperature TneLi (full triangles) and excited-state temper- 
atures (open symbols) as a function of the experimental excitation energy per nucleon 
(E )/(A Q ) for the reaction 197 Au on 197 Au at E/A = 1000 MeV. The indicated uncertain- 
ties are mainly of systematic origin. The displayed value of T^Be represents an average over 
the full range of (E )/(A ) (from Ref. H). 



quantum statistical model. The calculated TneLi exhibits a more continuous rise 
with decreasing than the thermodynamic temperature and is in very good 

agreement with the measured values. We thus find that, with the parameters needed 
to reproduce the observed charge partitions, this temperature-sensitive observable is 
well reproduced. A necessary requirement for a consistent statistical description of 
the spectator fragmentation is thus fulfilled. 



In the same experiment excited-state temperatures |3T], |32j were determined from 
the populations of particle-unstable resonances measured with the Si-Csl hodoscopes. 
The peak structures were identified by using the technique of correlation functions, 
and background corrections were based on results obtained for resonance-free pairs 
of fragments with Z < 3, such as p-d, d-d, up to 3 He- 7 Li. Correlated yields of p-t, 
p- 4 He, d- 3 He, 4 He- 4 He, and p- 7 Li coincidences and 4 He singles yields were used to 
deduce temperatures from the populations of states in 4 He (g.s.; group of three states 
at 20.21 MeV and higher), 5 Li (g.s.; 16.66 MeV), and 8 Be (3.04 MeV; group of five 
states at 17.64 MeV and higher). The probabilities for the coincident detection of 
the decay products of these resonances were calculated with a Monte-Carlo model 
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32|, |33|. The uncertainty of the background subtraction is the main contribution to 



the errors of the deduced temperatures. 

The values for the three excited-state temperatures are given in Fig. 8 as a 
function of the experimental excitation energy (Eq)/(Aq). Mutually consistent with 
each other, they appear to be virtually independent of the excitation energy, centering 
about a mean value of ~ 5 MeV. This is in striking contrast to the monotonically 
rising isotope temperature TheLi which is shown in comparison. 

A saturation of excited-state temperatures and a similar difference to the behavior 
of isotope temperatures has also been observed in central 197 Au + 197 Au collisions at 



incident energies E/A = 50 MeV to 200 MeV [34|. The interpretation given there 
starts from the fact that the excited states used for the temperature evaluation are 
very specific quantum states which may not exist in the nuclear medium in identical 
forms |35], |37j. The observed asymptotic states can develop and survive only 
at very low densities that may not be reached before the cluster is emitted into 
vacuum. Accordingly, the excited-state populations should reflect the temperature 
and its fluctuations at this final stage of fragment emission. The obtained mean value 
near 5 MeV is not inconsistent with results of dynamical calculations based on the 
BUU model [EDI. 



DENSITY 

Expansion is a basic conceptual feature of both the statistical multifragmentation 
and the liquid-gas phase transition. A volume of about six to eight times that oc- 
cupied at saturation density is assumed in the statistical multifragmentation models 
while the critical volume in the case of infinite nuclear matter has about three times 
the saturation value. The experimental confirmation of expansion or low breakup 
density is therefore of the highest significance. 

In central collisions of heavy nuclei, expansion is evident from the observation 



of radial collective flow fl38fl . Significant radial flow is not observed in spectator de- 
cays, and evidence for expansion has been obtained, indirectly, from model compar- 
isons. Models that assume sequential emission from the surfaces of nuclear systems 
at saturation density underpredict the fragment multiplicities while those assuming 
expanded breakup volumes yield satisfactory descriptions of the populated partition 



space [39]. The disappearance of the Coulomb peaks in the kinetic-energy spectra of 



emitted light particles and fragments, associated with increasing fragment produc- 
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Figure 9: Proton-proton correlation functions for the indicated intervals of Z^ound^ mea- 
sured at e iab 135° for the reaction 197 Au on 197 Au at E/A = 1000 MeV. The full lines 
are obtained from the analysis with the Koonin-Pratt formalism, the resulting source radii 
are indicated (from Ref. |lj]). 



tion, provides additional evidence consistent with volume emission or emission from 



expanded systems ED 



Interferometric methods permit experimental determinations of the breakup vol- 
ume or, more precisely, of the space-time extension of where the emitted products had 



suffered their last collision |yj. In the present case of spectator decay at relativistic 
energies the time scales should be rather short and the data, in good approximation, 
may be directly related to the breakup volume that is of interest here. 

Proton-proton correlation functions measured at angles of 0/ a & ~ 135° are shown 
in Fig. 9. They are characterized by a depression at small relative momenta, caused 
by Coulomb repulsion, and by a peak near q = 20 MeV/c that is caused by the 
S-wave nuclear interaction. Its comparatively small amplitude signals a large spatial 
extension of the proton source. The quantitative analysis of these data was per- 
formed with the Koonin-Pratt formalism f43fl . A uniform sphere was assumed for 
the proton source. The deduced radii are of the order of 8 to 9 fm which is dis- 
tinctly larger than the radius of 6.5 to 7 fm of a gold nucleus at normal density. 
The structure of the correlation functions is somewhat obscured at larger Z^ound but 
there is no indication that the radii and thus the volume should significantly change 
with impact parameter. The derived density, however, decreases considerably with 
increasing centrality, caused by the changing number of spectator constituents Aq 
(Fig. 10). These spectator masses result from the calorimetric analysis described 
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above and are found to be in good agreement with the prediction of the geometric 
participant-spectator model pi |. The mean relative density decreases to values be- 
low p/po = 0.2 for the most central bin, i.e. smallest Zb OU nd- The deduced values 
as well as the variation with centrality compare well with the densities entering the 
statistical multifragmentation model in the version that uses a fixed cracking dis- 
tance for the placement of fragments inside the breakup volume M. We recall here 
that the proton multiplicities and kinetic-energy spectra indicate a predominantly 
pre-breakup emission. This does not necessarily exclude the possibility that their 
interaction with the forming spectator matter causes the interferometric picture to 
reflect the extension of the latter. 

Besides the proton-proton correlations also correlations of other light charged 
particles were used to determine breakup radii. Pronounced resonances are exhibited 
by the p-a ( 5 Li, g.s.), d-a ( 6 Li, 2.19 MeV), and t-a ( 7 Li, 4.63 MeV) correlation 
functions. Their peak heights were analyzed using the numerical results of Boal and 



Shillcock [45|. The deduced values, within errors are in qualitative agreement with 
the proton values. Their much smaller error bars demonstrate that higher accuracies 
may be reached with these pronounced resonances of particle unbound states in 
light nuclei. The further development of the formalism needed for their quantitative 
interpretation seems therefore highly desirable. 



SUMMARY AND PERSPECTIVES 

New results for the mass, excitation energy, temperature, and density of excited 
spectator systems at breakup have been presented. The discussion of these data 
was meant to demonstrate that methods exist to determine these thermodynamic 
variables from the experiment. It was also intended to show that they are not without 
problems and that, perhaps, serious conceptual difficulties exist which may make it 
very difficult to arrive at unambiguous results. The problem caused for the excitation 
energy by the dynamical process of spectator formation is an example. 

It should be rewarding to further pursue this program, not only because of the 
hope to better identify signals of the liquid-gas phase transition but also because un- 
expected results may appear, as demonstrated for the temperatures. The saturation 
of the excited-state temperatures, according to the interpretation given here, is an 
interesting in-medium effect in itself. It also has the consequence that it may provide 
us with a means to determine rather reliably the internal temperatures of fragments 
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Figure 10: Spectator mass (Aq) (left panel) and relative density p/po (right panel) as a 
function of ^ound f° r the reaction 197 Au on 197 Au at E/A = 1000 MeV. The lines are 
meant to guide the eye. The densities were derived from proton-proton correlations (Fig. 
9) and from correlation functions constructed for light charged particles as indicated (open 
symbols, from Ref. [0]). 

at their final separation from the system. 

The interferometry with light charged particles confirms the low density of the 
breakup configuration. For more precise evaluations the formalisms needed to de- 
duce radii and densities from interferometric measurements will need continuing de- 
velopment. Measurements of radii for different fragment species, possibly emitted 
at different stages of the reaction, may allow us to test and refine the otherwise so 
successful picture of the single breakup state in the spectator decay. 
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